Strain sensitivity and superconducting properties of NbsSn from first principles calculations 
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Using calculations from first principles based on density functional theory we have studied the strain sensitiv- 
ity of the high-field superconducting magnet A15 Nb3Sn. The Nb3Sn lattice cell was deformed in the same way 
as observed experimentally on multi-filamentary, technological wires subject to loads applied along their axes. 
The phonon dispersion curves and electronic band structures along different high-symmetry directions in the 
Brillouin zone were calculated, at different levels of applied strain, e, both on the compressive and the tensile 
side. Starting from the calculated averaged phonon frequencies and electron-phonon coupling, the supercon- 
ducting characteristic critical temperature of the material, T c , has been calculated by means of the Allen-Dynes 
modification of the McMillan formula. As a result, the characteristic bell-shaped T c vs. £ curve, with a maxi- 
mum at zero intrinsic strain, and with a slight asymmetry between the tensile and compressive sides, has been 
obtained. These first-principle calculations thus show that the strain sensitivity of Nb3Sn has a microscopic 
and intrinsic origin, originating from shifts in the M>3Sn critical surface. In addition, our computations show 
that variations of superconducting properties of this compound are correlated to stress-induced changes in both 
the phononic and electronic properties. Finally, the strain function describing the strain sensitivity of Nb3Sn 
has been extracted from the computed 71(e) curve, and compared to experimental data from multi-filamentary, 
composite wires. Both curves show the expected bell-shaped behavior, but the strain sensitivity of the wire is 
enhanced with respect to the theoretical predictions of the bulk, perfectly binary and stoichiometric Nb3Sn. Un- 
derstanding the origin of this difference might open potential pathways towards the improvement of the strain 
tolerance in such systems. 

PACS numbers: 74.70.Ad, 74.25.Kc, 74.25.Jb, 68.35.Gy, 63.20.kd, 63.20.dk, 71.15.Mb 



I. INTRODUCTION 

The A15 phase NbsSn compound 1 is currently being used 
in a variety of large-scale scientific projects employing high- 
field superconducting magnets (above 10 T) 2 , including ITER 
(the International Thermonuclear Experimental Reactor) 3 " 5 , 
the 1 GHz NMR project 6 , and the CERN LHC Luminosity 
Upgrade 7 . In these high-field magnets, the mechanical loads 
during cooldown (due to different thermal contractions) and 
operation (due to Lorentz forces) can be very large, and since 
the superconducting properties of NbsSn strongly depend on 
strain 8 " 11 , an overall performance degradation can take place. 
Therefore, for a magnet's sound design it is of fundamental 
importance to have knowledge of the behavior of the super- 
conducting parameters (namely the critical temperature,T 
the upper critical field, B c 2, and the critical current, I c ) as a 
function of strain, e. 

Of particular interest is the uniaxial stress (either in ten- 
sion or compression), acting along the axial direction of the 
composite, multi-filamentary wires used in such systems: in 
a Cable-in-Conduit Conductor (CICC) 12 " 14 , for example, the 
NbsSn wires are inserted into stainless steel conduits, and 
compressive stresses due to the different thermal contraction 
coefficients of the different materials become important 15 ' 16 . 
Trasverse load components might also be important 11 ' 17 , but 
we will focus here only on the uniaxial ones. 

Within the framework of the Unified Scaling Law 18 , many 
authors 18 " 23 have proposed modified scaling equations which 
take into account the uniaxial strain dependence through the 



so-called strain function, s(e), but few attempts have been 
made 19 ' 20 ' 24 on obtaining a scaling law based on microscopic 
parameters. For this purpose, a very first step would be to ac- 
curately determine the electronic band structures, the phonon 
dispersion curves and the electron-phonon coupling terms, 
and study their evolution as a function of applied strain. To 
this aim, the knowledge of the NbsSn lattice cell deforma- 
tion when a multifilamentary wire is subject to different stress 




FIG. 1: The arrangement of atoms in the cubic (a) and tetragonal (b) 
phase of A15 Nb3Sn. The Nb atoms in the 2e sites form chains along 
the [001] direction, whereas those in the 4k sites are in chains along 
the [100] and [010] directions. To make the distortion clear, in (b), 
the cell is stretched to the abnormal value of £ = 40%, whereas the 
Poisson's ratio v is set to 0.4. 
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FIG. 2: The electronic band structure and density of states for M^Sn calculated for three representative strain states: +1.0% (•), zero applied 
strain (♦), and -1.0% The Fermi level is set to eV. 



components is of basic importance. Recent high resolution X- 
ray diffraction experiments on mechanically loaded samples 25 
have shown in detail how the Nb3Sn lattice cell deforms in 
the axial and the transverse directions; in particular, it was 
observed that the stress is completely transferred from the 
macroscopic level to the individual grains within the com- 
posite structure, so that a macroscopic uniaxial load directly 
corresponds to a stretching of the NbsSn lattice cell along 
the same direction, with the cell contracting in the transverse 
direction of an amount corresponding to a Poisson's ratio v 
equal to 0.38. 

The structural and electronic properties of Nb^Sn have been 
theoretically studied by several groups 26 " 31 , whereas the full 
phonon dispersion relations have been calculated by means of 
a tight-binding method 32 ' 33 and - more recently - by an ab ini- 
tio pseudopotential approach 34 . In particular, calculations by 
Tutiincii et al. give evidence of a strong interaction between 
the electronic states near the Fermi level and several phonon 
modes (longitudinal acoustic phonons and a group of optical 
phonon modes with average frequency of 4.5 THz) along the 
[111] direction. 

However, to the best of our knowledge no systematic ab ini- 
tio investigations have been made on studying the evolution of 
the band structure, phonon dispersion curves and supercon- 
ducting parameters (electron-phonon mass enhancement pa- 
rameter, X, and T e ) as a function of an applied uniaxial strain. 
In the present work, this issue has been addressed by em- 
ploying the plane-wave pseudopotential method, the density- 
functional theory, and a linear-response technique 35 ' 36 , and by 



using the results by Tutiincii et al. for the undistorted cell as a 
starting baseline for our calculations. 



II. DETAILS OF CALCULATIONS AND 
COMPUTATIONAL METHOD 

We used density functional theory and density-functional 
perturbation theory 36 as implemented in the QUANTUM- 
ESPRESSO software distribution 37 , within the local-density 
approximation 38 , a plane-wave expansion up to 40 Ry for the 
kinetic energy cutoff and ultrasoft pseudopotentials for Nb 
and Sn 39 . The Brillouin zone has been sampled on a 8x8x8 
Monkhorst-Pack (MP) mesh, corresponding to 126 special k- 
points within the irreducible part of the Brillouin zone (IBZ). 
We also have checked more dense grids (up to 16x16x16) but 
the results did not change considerably. 

Lattice dynamical calculations have been performed within 
the framework of the self-consistent density functional pertur- 
bation theory (DFPT) 36 , in which the dynamical matrices are 
calculated by sampling the IBZ with 8 independent q-points 
in the tetragonal phase. Dynamical matrices at any wave vec- 
tors can be Fourier deconvolved on this mesh, and the phonon 
dispersion curves along arbitrary symmetry directions can be 
easily obtained. In order to check the accuracy of the Fourier 
interpolation, we compared the results of this procedure with 
direct calculations on selected q-points not present in the grid. 

A denser grid of k-points (24x24x24 MP divisions) has 
been used in order to determine the electron-phonon interac- 
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FIG. 3: (a) Phonon dispersion curves of Nb3Sn at three different strain states: +1.0% (•), zero applied strain (♦), and -1.0% (»); (b) phonon 
DOS at each calculated strain. For completeness, experimental curves are also reported. 



tion parameter X, calculated as the Brillouin-zone average of 
the mode-resolved coupling strengths A. q; : 



(v = 0.4 25 ), and the distortions have been calculated in the 
strain range ± 1.0%, in steps of 0.2%. 



X=Y l W(q)X qj 
q; 



(1) 



where j indicates a phonon polarization branch, and W(q) 
are the weights associated with the phonon wavevectors q, 
normalized to 1 in the first Brillouin zone. 

The lattice parameter of the cubic cell is set to 5.29 A 40 ' 41 . 
For an uniaxial stress along the z-direction (a,) and under the 
assumption that the system is transversally isotropic, the strain 
state can be expressed as: 



e. 



: By ■■ 

E 



(2) 



which reflects the variation of the lattice parameters of the 
tetragonally distorted cell. In Eq. 2, v is the Poisson's ratio 
whereas E represents the Young's modulus. In our computa- 
tions, v has been set to the value measured in composite wires 



III. RESULTS AND DISCUSSION 

The Nb3Sn cubic phase belongs to the (P-i-3-) space 
group and the 0\ point-group symmetries, as shown in 
Fig. 1(a). The Sn atoms are situated on a bcc matrix whereas 
the faces of the cube are occupied by Nb atoms which form 
three sets of orthogonal chains along the principal axes. When 
a uniaxial strain is applied along the c-direction, the lattice is 
tetragonally distorted, with the Nb-chains in the [001] direc- 
tion differing from those in the [100] and [010] directions, as 
shown in Fig. 1(b). The distorted structure has a reduced 
symmetry Z40P^?). 

Starting from the cubic cell, a uniaxial strain has been ap- 
plied to the cell along the c-direction, according to Eq. (2). 
As a result, the deviatoric components of the strain lowered 
the system's point group symmetries: most of the phonon de- 
generacies have been removed and a changement in both the 
electronic and phonon dispersion bands have been induced. 
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FIG. 4: The behavior of: (a) the superconducting critical tempera- 
ture, T c ; (b) the el-ph coupling, X; and (c) the logarithmically av- 
eraged phonon frequency C0/n as a function of an applied uniaxial 
strain. Lines are guide for eye. 



A. Electronic band structures and phonon dispersion curves 

The calculated electronic structure along many high- 
symmetry directions of the simple-cubic Brillouin zone are 
displayed in Fig. 2 for three representative values of the ap- 
plied strain (zero, 1.0% and -1.0%). The energy bands of the 
cubic crystal are shown as diamonds, whereas the circles and 
the triangles represents the 1.0% and -1.0% strain states, re- 
spectively. Indeed, the tetragonal deformation does not affect 
the electronic structure in a severe way: the energy bands of 
the distorted and undistorted cell are almost unchanged, and 
there is no evident splitting of the cubic bands at the Fermi 
level, Ef. The electronic DOS is also not drastically affected 
by the tetragonal distortion. In Fig. 2, the Fermi level is 
marked by a dashed horizontal line and is set to eV. It is 
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FIG. 5: Density of states at Fermi level, N{Ef) calculated as a func- 
tion of strain. The curve is overlapped to the product function Xco 2 „, 
which in good approximation should be proportional to N(Ef) 42 . 



interesting to notice that Ep falls close to a sharp peak in the 
electronic DOS 27 , with a value for the density of states of the 
order of 20 states/eV. This peak is generated by several nearly 
dispersionless bands crossing the Fermi level in the F-M, 
T-R, and M-R directions and deriving from the 4d states of 
Nb atoms 34 . 

The phonon dispersion curves, calculated along several 
high-symmetry directions in the Brillouin zone, are plotted for 
three representative strains in Fig. 3; the phonon DOS is de- 
picted in the right panel. There is good agreement with previ- 
ous published results 34 at e = 0. For the sake of completeness, 
some inelastic neutron scattering measurements 43 " 45 are also 
reported in Fig. 3a, and again a good agreement is obtained. 



B. Derivation of superconducting T c as a function of strain 

The modifications of the phonon dispersion curves induced 
by a tetragonal distortion should have a strong effect on T c 
through the strain sensitivities of the averaged phonon fre- 
quencies C0/„, and of the el-ph coupling X. Therefore, both 
C0;„(e) and X(e) have been explicitly calculated, and r e (e) 
has been estimated by means of the Allen-Dynes modification 
of the McMillan formula 42 - 46 : 



hO)i„ '• 04 < 1+x > 

'JL e K-fi* (1+0.62*.) 

1.20 



(3) 



where /j* is the effective Coulomb-repulsion parameter 
which describes the interaction beetween electrons, and C0/„ 
is a weighted logaritmically averaged phonon frequency, de- 
fined as: 



(4) 
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where a 2 (co)F(co) is the Eliashberg spectral function. We 
assumed a negligible strain dependence of // compared to 
the other parameters, and frozen it as a constant in our 
computations 19 . Our results are reported in Fig. 4a-c. As far 
as the cubic phase is concerned, our finding are in good agree- 
ment with the ones by Tutiincu et al? A : a group of six phonon 
modes at the R-point (whose averaged phonon frequency is 
approximately 140 cm -1 ) are found to strongly interact with 
the p-d electronic states near the Fermi level. In these modes 
only the Nb chains vibrate, the Sn atoms being frozen at their 
equilibrium positions. The X^j corresponding to these modes 
are comprised in the range 0. 134-0.197. The overall electron- 
phonon interaction parameter {X = 1.85) agrees with the ex- 
perimentally measured value 47 , and - choosing /j* = 0.25 - the 
estimated critical temperature for the strain-free state is T c = 
18.3 K (also in agreement with the highest reported T c 48 ). 

As it can be clearly seen in Figs. 4b and 4c, both <B/„(e) 
and X(s) show a parabolic profile as a function of strain. Very 
close to the cubic phase (strain-free cell), C0/„ has a maxi- 
mum, implying a softening of the logarithmically averaged 
phonon frequencies when the system undergoes a distortion. 
The same behavior is found for X, whose maximum is ~ 1.85. 
The strength of the el-ph interaction weakens as |e| increases. 
As a result, by varying the axial strain the T c curve assumes the 
characteristic bell shape (Fig. 4a) 49,50 . In addition, the curve 
shows a slight asymmetry with respect to the maximum, due 
mainly to an asymmetry in the phononic contribution (C0/„). 
However, a clear confirmation of this would deserve a more 
detailed analysis, based on an increased density of points on 
the curve. Qualitatively, the T c vs. e curve reproduces the 
experimental strain sensitivity of the critical current found in 
A15 superconductors: as it is well known, when a Nb3Sn mul- 
tifilamentary wire is subjected to a longitudinal strain, its crit- 
ical current shows a maximum at zero intrinsic strain, and de- 
creases reversibly with the applied load (see for example Ref. 
51 and references therein), with a slight asymmetry 52 between 
the compressive and tensile sides. In this sense, these first- 
principle calculations suggest that the origin of such strain 
sensitivity in Nb3Sn is intrinsic and microscopic in its nature. 

Our calculations also show that N(Ep) is influenced by 
strain. This quantity is related to the el-ph coupling constant 
and to the averaged phonon frequency trhough the following 
expression 42 : 

N(E F )=M<I 2 >(x) 2 RMS X (5) 

in which < I 2 > is the average over the Fermi surface of 
the el-ph matrix element squared, (x)rms is a weighted RMS 
phonon frequency and M is the average ionic mass. By fur- 
ther assuming that the strain sensitivity of the normalized av- 
eraged frequencies Wrms an d <*>/„ are the same 54 and that 
< I 2 > does not depend on any applied strain, it follows that 
N(Ef,e) oc X(ti 2 n . The strain dependence of N(Ef), calcu- 
lated either directly or through Eq. 5 are consistent with one 
another, as can be clearly seen in Fig. 5. Results of Figs .4 
and 5 thus provide evidence that the strain is affecting both 
the phononic and the electronic properties in the same way. 
Indeed, existing studies consider only the two extreme cases 
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FIG. 6: Direct comparison between the theoretical (dashed line) and 
experimental (line and markers) strain function 53 . The theoretical 
s(e) correctly reproduces a bell shape, the mismatch being attributed 
to those extrinsic effects generally observed in technological wires. 



where either the lattice deformations or the electronic prop- 
erties modifications are considered as a source for the strain 
sensitivity. Markiewicz 23 ' 24 ' 55 has attempted to correlate the 
microscopic full invariant analysis with the Elisahberg-based 
relations for T c through strain-induced modifications in the 
electron-phonon spectrum. In this model, however, the strain 
induced changes in N(Ep) are only accounted for through 
a strain-modified frequency dependence of the el-ph interac- 
tion. In other words, the changes in N(Ep) are not directly 
calculated, although the model is sufficiently accurate. Fur- 
thermore, microscopical theoretical predictions by Taylor and 
Hampshire 19 and Markiewicz 55 have shown that the variation 
of the superconducting properties of Nb3Sn multifilamentary 
wires submitted to uniaxial strain are correlated to changes of 
the phonon spectrum, rather than to the electronic density of 
states. On the other side, many works 54 ' 56 ' 57 have linked the 
superconducting properties of A15 compounds to the varia- 
tions in the electronic properties and N(Ef). 

From the experimental point of view, it is important to 
highlight new methods for extracting the electron DOS from 
resistivity data have been explored with the aim to study 
whether the strain sensitivity is correlated to the electronic 
modifications 58 . 

However, according to our results, any model aiming to 
describe the superconducting properties of A15 compounds 
from microscopic theories should take both contributions into 
account. 



C. Comparison with experimental results 

Starting from the pioneering work of Ekin 18 , in many mod- 
els available in literature 18 " 23 the strain sensitivity of Nb3Sn 
is generally parameterized using the strain function, s(e), de- 
fined as 18 ' 59 : 
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S (e) = 



B c2 (e,T = 0) 
B c2 (0,T = Q) 



Tc(0)\ 



(6) 



where B C 2(£,T) represents the superconducting upper crit- 
ical field, depending on strain and temperature, and where w 
f» 3 for A15 materials 18 . 

The strain function calculated using Eq. 6 is plotted in 
Fig. 6 (dashed line) and compared with the curve extracted 
from experimental I c vs.e measurements 53 on a technological, 
multifilamentary wire (lines and markers in the figure). Al- 
though both curves show the expected bell shape, the experi- 
mental and theoretical curves are quite different. In particular, 
the strain sensitivity of the composite wire is enhanced when 
compared to the theoretical prediction of a bulk, perfectly bi- 
nary and stoichiometric Nb3Sn system, and this remains valid 
also if the value of w is increased (within physical accepted 
ranges). This is not surprising, and the difference might 
have either an intrinsic or extrinsic origin. Among extrinsic 
phenomena inducing performance degradation with strain in 
technological wires 60 , filament breakage, reduction in wire's 
cross-sectional area, stress-induced martensitic formation 53 , 
and microcrack/defect formation in the superconductor might 
play a role. The reversibility of the experimental data plotted 
in Fig. 6 implies that some of these phenomena can be ruled 
out, e.g. filament breakage, microcrack and extended defect 
formation. All others mentioned above, being reversible, can 
in principle explain the strain sensitivity in NbsSn, but their 
effect over I c is small and cannot account for the observed 
behavior 60 . As far as intrinsic mechanisms are considered, 
our calculations unambiguously show that strain sensitivity 
in Nb3Sn is associated to lattice and electronic deformations, 
which result in shifts in the M^Sn critical surface. 

It should be underlined the fact that in this first principles 
study we have neglected any sublattice displacements of the 
Nb atoms 27 leading to Nb-chains dimerization. This can also 
have an effect on the theoretical s(e). In our calculations, Nb 
atoms are frozen in their ideal positions, and the cubic cell is 



stable with respect to a tetragonal strain 27 . However, if the 
Nb atoms are allowed to relax, a dimerization of the chains 
occurs, and the undistorted cubic structure become unstable 
with respect to a spontaneous sublattice distortion. However, 
such distortions are small (the ratio between the major and mi- 
nor axis of the tetragonal cell spans from 0.9938 to 0.9964 27 ) 
and therefore their effect on s(e) is expected to be negligible. 
Also, our system is perfectly binary, whereas in technological 
wires Ti or Ta is inserted as the ternary element, with the aim 
to improve the pinning efficiency and therefore J C (B). More- 
over, due to compositional inhomogeneities, a distribution of 
the superconducting properties is observed (T c depends on the 
atomic Sn content, having a maximum at Sn% = 0.25%); see 
for example Refs. 61 and 62 and reference therein. Consider- 
ing all these aspects, it is clear that the polycrystalline M^Sn 
formed by a reaction heat treatment inside a composite, multi- 
filamentary system has a microscopic structure that inevitably 
deviates from that of an ideal NbsSn lattice cell. Therefore, 
differences between the theoretical computations of s(e) and 
the experimental degradation in wires can be expected, which 
might possibly suggest paths towards the improvement of the 
strain tolerance in technological wires. For example, at 0.5% 
compressive strain the strain sensitivity could in principle be 
reduced to ~ 0.07%, thus helping in the design of those de- 
vices where the levels of strain experienced by Nb3Sn are 
sufficiently large, as is the case of high-field superconducting 
magnets. 



Acknowledgments 

We thank Federico Quagliata and Pietro D'Angelo for set- 
ting up the environment and the CRESCO parallel cluster at 
ENEA C.R. Frascati. M.B.N, wishes to acknowledge partial 
support from the Office of Basic Energy Sciences, U.S. De- 
partment of Energy at Oak Ridge National Laboratory under 
Contract No. DE-AC05-00OR22725 with UT-Battelle, LLC. 



* demarzi@enea.it 

1 B. T. Matthias, Phys. Rev. 92, 874 (1953). 

2 T. Miyazaki, T. Hase, and T. Miyatake, Handbook of Supercon- 
ducting Materials (Bristol: IOP Publishing, 2003), pp. 639-72. 

3 A. Vostner and E. Salpietro, Supercond. Sci. Technol. 19, S90 
(2006). 

4 N. Mitchell, A. Devred, P. Libeyre, B. Lim, S. E, and the ITER 
Magnet Division, IEEE Trans, on Appl. Supercond. 22, 4200809 
(2012). 

5 A. Devred, I. Backbier, D. Bessette, G. Bevillard, M. Gardner, 
M. Jewell, N. Mitchell, I. Pong, and A. Vostner, IEEE Trans, on 
Appl. Supercond. 22, 4804909 (2012). 

6 H. Wada and T. Kiyoshi, IEEE Trans. Appl. Supercond. 12, 715 
(2002). 

7 L. Bottura, G. de Rijk, L. Rossi, and E. Todesco, IEEE Trans, on 
Appl. Supercond. 22, 4002008 (2012). 



s G. Rupp, I. Appl. Phys. 48, 3858 (1977). 
9 J. Ekin, IEEE Trans. Magn. 15, 197 (1979). 

10 B. ten Haken, Ph.D. Thesis, Technical University of Twente, En- 
schede, The Netherlands (1994). 

11 J. Ekin, J. Appl. Phys. 62, 4829 (1987). 

12 M. O. Hoenig, A. G. Montgomery, and S. J. Waldman, IEEE 
Trans, on Mag. 15, 792 (1979). 

13 P. Bruzzone, IEEE Trans. Appl. Supercond. 16, 839 (2006). 

14 M. Spadoni, IEEE Trans, on Magn. 30, 1699 (1994). 

15 A. Nijhuis and Y. Ilyin, Supercond. Sci. Technol. 19, 945 (2006). 

16 N. Mitchell, Supercond. Sci. Technol. 18, S396 (2005). 

1 7 

G. Mondonico, B. Seeber, A. Ferreira, B. Bordini, L. Oberli, 
L. Bottura, A. Ballarino, R. Fliikiger, and C. Senatore, Supercond. 
Sci. Technol. 25, 115002 (2012). 

18 J. W. Ekin, Cryogenics 20, 61 1 (1980). 

19 D. M. Taylor and D. P. Hampshire, Supercond. Sci. Technol. 18, 



7 



40 



S241 (2005). 

S. Oh and K. Kim, J. Appl. Phys. 99, 0330909 (2006). 

A. Godeke, B. ten Katen, H. H. J. ten Kate, and D. C. Larbalestier, 
Supercond. Sci. Technol. 19, R100 (2006). 

D. Arbelaez, A. Godeke, and S. Prestemon, Supercond. Sci. Tech- 
nol. 22 (2009). 

W. D. Markiewicz, Cryogenics 46, 864 (2006). 

W. D. Markiewicz, Cryogenics 44, 767 (2004). 

L. Muzzi, V. Corato, A. della Corte, G. D. Marzi, T. Spina, 

J. Daniels, M. D. Michiel, F. Buta, G. Mondonico, B. Seeber, 

et al., Supercond. Sci. Technol. 25, 054006 (2012). 

B. M. Klein and Z. Lu, Physica B: Condensed Matter 296, 120 
(2001). 

B. Sadigh and V. Ozolins, Phys. Rev. B 57, 2793 (1998). 

Z. W. Lu and B. M. Klein, Phys. Rev. Lett. 79, 1361 (1997). 

L. F. Mattheiss and W. Weber, Phys. Rev. B 25, 2248 (1982). 

B. M. Klein, L. L. Boyer, D. A. Papaconstantopoulos, and L. F. 

Mattheiss, Phys. Rev. B 18, 641 1 (1978). 

L. F. Mattheiss, Phys. Rev. B 12, 2161 (1975). 

W. Weber, Physica B+C 126, 217 (1984), ISSN 0378-4363. 

W. Weber, Electronic Structure of Complex Systems (Plenum 

Press, New York, 1984), vol. 113, p. 345. 

H. M. Tiitiincii, G. P. Srivastava, S. Bagc l, and S. Duman, Phys. 
Rev. B 74, 212506 (2006). 

S. Baroni, P. Giannozzi, and A. Testa, Phys. Rev. Lett. 58, 1861 
(1987). 

S. Baroni, S. de Gironcoli, A. Dal Corso, and P. Giannozzi, Rev. 
Mod. Phys. 73, 515 (2001). 

QUANTUM-ESPRESSO is a community project for high-quality 
quantum-simulation software, based on density-functional theory, 
and coordinated by Paolo Giannozzi. See http://www.quantum- 
espresso.org and http://www.pwscf.org. 
N. Troullier and J. L. Martins, Phys. Rev. B 43, 1993 (1991). 
We used the pseudopotentials Nb.pw91-nsp-van.UPF and 
Sn.pw91-n-van.UPF from the http://www.quantum-espresso.org 
distribution. 

R. G. Maier, Z. Naturforsch. Teil A 24, 1033 (1969). 



41 

42 
43 

44 
45 
46 
47 

48 
49 

50 

51 

52 
53 

54 

55 
56 
57 

58 

59 
60 
61 

62 



V. Guritanu, W. Goldacker, F. Bouquet, Y. Wang, R. Lortz, 

G. Goll, and A. Junod, Phys. Rev. B 70, 184526 (2004). 

W. L. McMillan, Phys. Rev. 167, 331 (1968). 

L. Pintschovius, H. Takei, and N. Toyota, Phys. Rev. Lett. 54, 

1260 (1985). 

J. D. Axe and G. Shirane, Phys. Rev. B 8, 1965 (1973). 
G. Shirane and J. D. Axe, Phys. Rev. B 18, 3742 (1978). 
P. B. Allen and R. C. Dynes, Phys. Rev. B 12, 905 (1975). 
E. L. Wolf, J. Zasadzinski, G. B. Arnold, D. F. Moore, J. M. Row- 
ell, and M. R. Beasley, Phys. Rev. B 22, 1214 (1980). 
J. Hanak, K. Strater, and R. Cullen, RCA Review 25, 342 (1964). 
R. Fliikiger, R. Isernhage, W. Goldhacker, and W. Specking, Adv. 
Cryo. Eng. (Materials) 30, 851 (1984). 

J. W. Ekin, Experimental Techniques for Low -Temperature Mea- 
surements (New York: Oxford University Press, 2007). 

A. Godeke, Ph.D. Thesis, Technical University of Twente, En- 
schede, The Netherlands (2005). 

R. Fliikiger, D. Uglietti, V. Abacherli, and B. Seeber, Supercond. 
Sci. Technol. 18, S416 (2005). 

G. D. Marzi, V. Corato, L. Muzzi, A. della Corte, G. Mondonico, 

B. Seeber, and C. Senatore, Supercond. Sci. Technol. 25, 025015 
(2012). 

K. C. Lim, J. D. Thompson, and G. W. Webb, Phys. Rev. B 27, 
2781 (1983). 

W. D. Markiewicz, Cryogenics 44, 895 (2004). 
W. Weber and I. Goldberg (????). 

B. M. Klein, L. L. Boyer, D. A. Papaconstantopoulos, and L. F. 
Mattheiss, Phys. Rev. B 18, 641 1 (1978). 

M. G. T. Mentink, M. M. J. Dhalle, D. R. Dietderich, A. Godeke, 

W. Goldacker, F. Hellman, and H. H. J. ten Kate, AIP Conference 

Proceedings 1435, 225 (2012). 

D. O. Welch, Adv. Cryo. Eng. 26, 48 (1980). 

J. Ekin, IEEE Trans. Magn. 13, 127 (1977). 

C. Senatore, V. Abacherli, M. Cantoni, and R. Fliikiger, Super- 
cond. Sci. Technol. 20, S217 (2007). 

A. Godeke, Supercond. Sci. Technol. 19, R68 (2006). 



